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Abstract 

Ionic currents accompanying DNA translocation strongly depend on molarity of the electrolyte 
solution and the shape and surface charge of the nanopore. By means of the Poisson-Nernst-Planck 
equations it is shown how conductance is modulated by the presence of the DNA intruder and as 
a result of competing electrostatic and confinement factors. The theoretical results reproduce 
quantitatively the experimental ones and are summarized in a conductance diagram that allows 
distinguishing the region of reduced conductivity from the region of enhanced conductivity as a 
function of molarity and the pore dimension. 
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I. INTRODUCTION 


Translocation of DNA in narrow pores plays a central role in many biological processes, 
snch as viral infection by phages and interbacterial DNA transdnction pQ as mnch as in the 
development of novel devices for high-thronghpnt and low-cost biotechnological applications 
In the last decade, several experimental studies have explored the translocation process 
of DNA through protein channels across cellular membranes or microfabricated channels |1] . 

Solid state nanopores have attracted considerable attention as ideal devices for reading 
off the base pairs while tracking the DNA motion through nanopores [S]. Experimentally, 
during the DNA translocation through a nanopore one measures rapid variations of the ionic 
currents due both to the steric hindrance and the negative DNA charge density distribution. 
In principle, these variations can be used to decypher the DNA base pair sequence, if one 
is able to finely identify the factors which modulate the current. To this purpose, one 
must characterize the differences between the ionic signal when the DNA is in the pore or 
away from it. One observes that the ionic current shows a larger/lower value when the 
DNA is inside/outside the pore depending on the concentration of the electrolyte. Since the 
reservoirs provide an amount of counterions sufficient to screen the DNA charge the number 
of charge carriers increases and thus determines a larger conductance, in spite of the fact 
that the effective section available to the passage of ions is smaller. 

From the theoretical point of view DNA translocation and the accompanying electroki- 
netic transport of electrolytes involve the comprehension of the competition between elec¬ 
trostatic, excluded volume and fluid-atom hydrodynamic couplings, with the key role played 
by the local confinement. Therefore, understanding the physical mechanisms that regulate 
ionic transport calls for an accurate determination of each contribution and, wherever pos¬ 
sible, for realistic computational modeling. In addition, as the diameters of the pore and of 
DNA can be as small as a few nanometers, ionic transport is genuinely microscopic. At such 
length scale a careful evaluation of the local interactions requires a sophisticated theoretical 
treatment [6]. 

When DNA translocates in a solid state pore, it assumes an elongated conformation, due 
to both the electrostatic repulsion stemming from different parts of the DNA backbone, 
and by the axisymmetric and narrow shape of the pore. One of the most interesting effects 
of the high confinement is that the ensuing ionic current can be modulated in different 
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and opposing ways. In fact, several anthors have reported that ionic cnrrent can be either 
enhanced or blocked by the presence of DNA within the pore [M2]. Cnrrent blockage refers 
to the fact that ionic cnrrent conld be temporarily rednced when DNA is present within the 
channel. Snch effect is typically associated with the occlusion due to DNA that diminishes 
the ionic fluxes. Vice versa, current enhancement is also possible in a low-concentration 
electrolyte. Such effect can be ascribed to the excess of charge carriers accompanying DNA 
as it occupies the pore region, that are loosely bound to the biomolecule and therefore 
available to conduction. Both blockage and enhancement are observable in experiments and 
their occurrence depends on the molarity, the degree of conhnement, the nanopore material 
(such as SiN or Si02), and possibly other physico-chemical parameters [2l 113]. 

In the present paper, we analyze how current modulation depends by several parameters 
by employing a dual theoretical/computational approach. Due to the high level of conhne¬ 
ment, one usually ignores the presence of convective currents, that is, considers frictional 
forces strong enough to effectively damp out electro-osmotic effects. To test such an issue 
we shall compare results obtained by taking into account convective currents with those 
that exclude these contributions. Our description is encoded by the well-known Poisson- 
Nernst-Planck (PNP) theory [T3], as largely employed in the study of ionic transport in 
biological ion channels, such that only electrostatic, diffusive and migration contributions 
are considered. In addition, in order to simplify the analysis an effective one-dimensional 
equation based on a suitable homogenization procedure is written down [T3] . By identifying 
the entropic forces arising from the degree of conhnement and, by imposing the local elec¬ 
troneutrality condition, the conductance of the pore-DNA system can be estimated. The 
enhancement/blockage diagram of the system reveals how the ionic current is modulated as 
a function of molarity, pore geometry, ehective charge on DNA and pore surface. In addi¬ 
tion, besides the role of the DNA intruder in the pore, the equations allow determining the 
contribution of the reservoirs, with the related access resistance, to the overall transport. 
Being a one-dimensional diherential equation, the analysis is particularly manageable in 
terms of computational and modeling ehorts. For instance, we will represent DNA as either 
a smooth or a corrugated cylinder, being either at rest or in motion within the pore, and 
assess how the enhancement/blockage diagram varies correspondingly. With respect to our 
previous study [13], we have improved the solution by replacing the local electroneutrality 
approximation used to determine the electric field by the more rehned solution of the Pois- 
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son equation. The computational cost of the approach is modest and can serve as a useful 
method to perform an early survey of DNA/pore systems. 

The paper is organized as follows in section we dehne the model and its governing one 
dimensional effective equations. In section [ITT] we perform numerical calculations and discuss 
the results concerning the conductance and a conductance diagram obtained by considering 
the variation of the conductance upon inserting into the pore a long thin cylindrical charged 


object mimicking the DNA molecule. Finally, in section W we present a few concluding 
remarks. 


II. MODEL SYSTEM AND REDUCTION TO A ONE-DIMENSIONAL EFFEC¬ 
TIVE PROBLEM 



X — Iq X = l\ X = I2 


X —I3 X — U X — Is 


FIG. 1: Geometry of the pore region. 

We start by describing the system and how it is modelled in effective one-dimensional terms. 
The geometry of the three dimensional pore is sketched in Fig. and consists of two 
conical funnels connecting a straight pore. The radius of the channel, which has cylindrical 
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symmetry, can be written as a function of x alone: 


Rrh{x) = 


-Ro 

Ro - Ro(a; 
Ri 

Ro + Bq[x 


li Iq < X < li 
h) h < X < I2 
ii I2 < X < 

I4) ii I3 < X < I4 


( 1 ) 


Ro if ^4 < a; < I 5 

\ 

The two conical funnels of variable radius are joined by the channel with section of constant 
radius Ri = Rq — Ro (^2 ~ ^ 1 ) and length L = I 3 — I 2 , whose inner walls carry a hxed surface 
charge of area density Eg, Rq is the slope of the funnel and moreover d = I 2 — I 1 = ^4 — /s and 
6 = h —If) = I 5 — 14 . Each funnel connects the pore to a reservoir and / 2,3 are the coordinates 
of the inlet and outlet of the cylinder, respectively. The DNA molecule is described as a 
very long thin cylinder of radius Rdna and with Rdna{x) < Ri < Rq, whose axis coincides 
with the one of the pore [15]. The space between the pore walls and the DNA is hlled by 
an electrolytic solution and the cathode and the anode are placed at x = Iq and x = I 5 , 
respectively, and a potential difference is applied between them, as shown in Fig. 

The evolution of the ionic concentrations is described within the framework of the Poisson- 
Nernst-Planck theory, which combines Pick’s diffusion law with the drift induced by the 
electric held and the convection of the ions due to the motion of the solvent determined 
by the DNA movement [12]. The PNP three-dimensional model can be solved numerically, 
however an effective one-dimensional description can shed some light and has gained increas¬ 
ing popularity in recent times [miiE]. In situations of interest the geometry of the channels 
is such that it is possible to neglect the variations of the relevant helds along some direc¬ 
tions and reduce the description to a one-dimensional problem. In nanochannels of constant 
section, whose transverse width is much shorter than the longitudinal size, the reduction is 
straightforward, otherwise one must take into account the effect of the varying shape. 

The lengths which determine the electrostatic properties of the system are the Debye 
screening length , Xd, the Bjierrum length Xb, the typical radius of the pore, of DNA and the 
pore length. The Debye length is a function of the bulk densities of each ionic species, Ub, the 
dielectric permittivity, e, and the temperature T, the electronic valence z according to Ad = 
where ks is the Boltzmann constant and e the elctronic charge. The Bjerrum 
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length is Xb = under standard ambient conditions and for aqueous solutions is 

typically of the order of 0.7 nm, while Xb in electrolytic solutions of concentrations in the 
range 0.1 — 1.0 M, Xb varies between ~ 0.961 and 0.304 nm. Thus in channels of nanometric 
diameter Xb can exceed their transverse size and the ensuing behaviour of the solution is 
mainly determined by the presence of the surface charges due to the partial overlap of the 
double layers of different bounding surfaces. One hnds that the resulting ionic atmosphere 
is mainly composed by counterions, a feature of capital importance because it yields the 
possibility of controlling through the surface charge the ionic currents through the pores. 

In the present treatment we only consider the steady current regime which does not re¬ 
quire the full time dependent solution of the PNP equations. The condition that the ionic 
current densities, J^, must have vanishing divergence, V ■ J^(r) = 0, together with the 
property of impenetrability of the walls, yield the following relation between the average 
axial component, Jx, of the three-dimensional current and the variable section of the pore: 

= {J^{x))S{x) = constant (2) 


where *S'(x) is the transverse section of the system and the current. After this premise, 
in the stationary state /“ is the sum of four different terms: the diffusive contribution, 
the entropic term, accounting for the modulation of the conhnement [I9H2I], the migration 
contribution stemming from the driving electric held and a convective term due to the 
electro-osmotic how induced by the charges present on the pore and DNA surfaces. The 
r.h.s. of the following expression encodes the four diherent contributions. 


(ic^(x) dlnS{x) 


dx 


dx 


N ez^{E{x)) ^ 
c^{x) -—— c^{x) 


-|-c^(x) 


ksT 

'^conv 




(3) 


where is the species coefficient of dihusion. 

The convective velocity Vconv, will be discussed below, depends on the surface charges 
of the pore and DNA and on the electric held along the axis direction. Equation ([^ 
represents an ordinary hrst order diherential equation for the one-dimensional linear density 
of species ±, which is related to the three dimensional ionic densities n^{r,t) by the 
following sectional averaging: 


c^(x, f) = (n^(x, t))S'(x) = f dSn^{r,t) 

Js{x) 


(4) 
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where the sectional area S{x) for the system shown in Fig. [^is 

S{x) = Tr(^Rlf,{x) - Rdna{x)) . 
Similarly, we dehne the sectionally averaged electric field: 

1 i 


{E{xR)) = 


dSE^{r,t ), 


J S{x) 

which as shown in ref. [I51E2] satishes the following differential equation 

d{E{x)) ^ dlnS{x) 

1 


( 5 ) 


( 6 ) 


■{E{x)) 


eS{x) 


dx dx 

z'^c^{x) + z~c~(x) + 27iRch{x)Tiu,{x) 


X 


l+l^!^y + 2,rRo^^EoMA 


(7) 


where S^(a;) is the wall charge density, being Sq in the pore and zero elsewhere. 

The terms proportional to the surface charge of the walls Sq and of the DNA intruder 
Edna appear as source terms in the one-dimensional representation and stem from the 
boundary conditions and on the reduction of the original three dimensional problem. The 
geometrical term containing the logarithm of S{x) enforces the conservation of the flux of 
the electric held in a system of variable pore section. 

The set of effective one-dimensional equations here considered represents a convenient 
approach to the solution of full three-dimensional problems, which involve a demanding 
numerical effort. As far as the performance of the one-dimensional PNP model (PNPl) 
versus the corresponding three dimensional model (PNP3) we verihed that the narrower the 
channel, the larger is the difference in prohles between the two types of solutions [6]. 

In the case of a binary electrolytic solution one has to solve three coupled differential 
equations: the two equations ([^ are needed to determine the local concentrations of the 
± species, while the third equation ([^ determines the electric field produced by the ionic 
charge distribution and the hxed charges. In the literature one often employs the so-called 
local electroneutrality approximation (LEN) |2S1123 because it avoids the effort of Ending 
the solution of eq. ([^ which is not known in closed analytical form. The LEN approximation 
states that locally there is no charge separation, a fact that may be justified because the 
typical values of the screening length Xd are small as compared with the system dimensions, 
such as the longitudinal and transverse channel sizes. In other words one assumes that at 
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each point x the mobile charges exactly balance the fixed surface charges, as it occurs in 
the limit of vanishing A/j- The electric field does not satisfy the one dimensional Poisson 
equation (j^, but takes values that enforce the absence of charge separation. Under steady 
state conditions one can express the charge density distribution in terms of the fixed surface 
charge distribution and obtain a very simple equation for the density profiles in terms of the 
electric current and the mass current ng. On the other hand, dimensional analysis shows 
that the LEN may break down in nanometric systems, that is, when the system size is of 
the same order as 

To fix the convective velocity, Vconv, that plays the role of an input parameter in the 
PNP model, we follow Ghosal’s treatment [25H2Z]. The presence of surface charges and of 
the applied electric field, E, along the axis, generates an electroosmotic flow (EOF), whose 
velocity is parallel to E. The resulting radial velocity profile, u{r) of the fluid, vanishing at 
the pore surface, reads: 

Uq (j)w 

where the potential 0(r) is the solution of the two dimensional Debye-Huckel equation in 
the anular region Rdna < ^ < Ri- 


0(r) = 


^DNA^D 


D 


AIo —]+BKo — 


A 


D 


(9) 


(f)w is the value at r = Ri, uq = eEcl)^,/fi a characteristic electroosmotic velocity and p the 
dynamic viscosity of the solution. Kq and Jq are modified Bessel functions of integral order. 
The constants A and B are determined (see ref. [25]) by imposing the boundary conditions 


— e4>'{RDNA) = '^DNA ( 10 ) 

e4>'{Ri) = So. (11) 

We approximate the convective velocity featuring in eq. ([^ by the following sectional 
average: 

2 

Vconv = p 2 _ p 2 - / drru{r). (12) 

■^1 ^DNA J Rdna 

Notice that in the cases here considered such a velocity turns out to be larger than the 
translocation velocity of the DNA molecule whose value can be estimated as: 

4>{Rdna) — 4>w 

VdNA — Uo - 7 -. 

4^w 
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In order to give an idea of the importance of the electric current associated with the EOF, 
I EOF-, with respect to the conduction current, lohm, we consider the ratio : 

It- 


leof _ ksT 2x2 X 
- - 


lohm ^ jdD 

For sodium ions in water at 1 M concentration the above ratio is approximately 0.003, 
where D and p are the diffusion coefficient of the ions and the dynamic viscosity of water, 
respectively. 

As mentioned in the introduction to detect the passage and eventually sequence DNA, one 
monitors the ionic conductance G of the pore. On theoretical grounds, the total resistance 
^tot of the system shown in Fig. 1 is the sum of the resistances of five different pieces: 
the resistance associated with the two funnels plus the resistance due to the two cylindrical 
regions adjacent the electrodes and the central part, ^tot = ^ch + 23fJcone + “I^reservoir- In 
the absence of surface charges the total conductance can be approximated by the following 
Ohmic formula (see ref. [TJ]) : 

71 1 


G = 


n L 2S 

mi ' Bo^Ri Ro- 


^ Bn \ R, Rn 7 


kaT 1 


(13) 


-. We first 


where po is the resistivity defined in terms of the coefficient D as po = 
evaluate the effect of the presence of a cylindrical intruder of radius Rdna within the pore 
in the limit of i?o >> Ri nnd L » Ri and no charges. By taking into account the 
reduced section available one estimates the difference in conductance between the free and 
the obstructed pore we find from eq. ( [I^ : 

71 1 


G{R{)-G{^1 


2 _ f?2 \ ~ 

^DNA I — ^ r ^DNA■ 
Po E 


(14) 


Although the conductance due to the obstruction always results negative on the basis of steric 


arguments, the presence of negative charges on the DNA and not included in formula (13), 


may lead to a different answer. In fact, some positive ions migrate from the reservoirs towards 
the negative surfaces and increase locally the number of mobile carriers and one observes a 
conductivity larger than in the absence of DNA. Such a phenomenon is particularly relevant 
at low ionic concentrations and is related to the so-called surface conduction mechanism 


The competition between enhanced conduction due to the presence of surface charges 
and the depletion due to the reduction of the available pore section due to the DNA leads 
to an interesting dependence of the relative conductance on the DNA features. 
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Hereafter, we shall study the conductance variations, as the concentration of the solution 
and the surface charge are varied, by comparing the numerical results relative to the one¬ 
dimensional PNP model with the predictions based on the LEN theory, which provides 
simple expressions for small potential and concentration drops. 

III. METHODS AND NUMERICAL RESULTS 

The simulations were performed considering systems comprised of the hve different regions 
illustrated in Fig. [Tj In the numerical calculations Ri took values between 5nm and 
20 nm, Bq = (i?o ~ Ri)/d while the other geometrical parameters were always Rq = 25nm 
L = SAnm, 6 = 3nm , d = 20 nm. The DNA is modelled as a long rigid cylinder of radius 
Ron A = 1.1 um, coaxial to the nanochannel, and having surface charge density 
This schematic representation of DNA is justihed by the persistence length of the dsDNA 
molecule which is about 50 nm. The surface charge densities were assumed to have values 
So = —0.375 e/nm^, S^wa = —0.38175 e/nm^, Rdna = 1.1 nm. A voltage drop was applied 
on the ends of the system: 

0(/o) = (t)i = (t)o- (15) 

In order to solve numerically eqs. (|^ and Q we have introduced a one-dimensional mesh 
Ax and dehned non dimensional quantities in the following way: the Debye length, k~^ = 
k~^Ax^ the concentration, n = nAx^, the charge density S = SAa:^/e, the electric potential 
■0 = e'ljj/kBT, the applied potential difference AV = eAU/fcsT. 

In the following, unless explicitly stated we shall perform calculations by neglecting the 
convective contribution to the current and include this term only when computing the global 
phase diagram. 

We started by measuring the dependence of the total electric current. Rot = I~ — , on 

the applied potential difference AV for various values of the salt concentration relative to a 
solution of ions of identical masses and mobilities and extracted the conductance G. In the 
case of the empty pore we found the I—V characteristic to be linear at concentrations 0.01 M, 
0.1M,1.0M (data not reported). On the other hand, the conductance as concentration is 
lowered at hrst decreases almost linearly, showing a standard bulk Drude behavior, but 
hnally reaches a plateau value when is of order of the ratio So/7?i. 

In order to understand the effect of the two vestibules we compared the results relative 
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to the set up of Fig. [TJwith those relative to a system with constant radius Rq = Ri, that 
is a straight cylindrical pore. Fig. [^displays the ionic current versus molarity in these two 
cases for the empty system and for the system partially occluded by DNA when Ri = 5 nm. 
We observe that the total ionic current in the non uniform channel, by virtue of its larger 
cross section, is larger than the one relative to the straight model. In the inset of Fig. 
the currents of counterions and coions are reported separately, limited to the case of the 
inhomogeneous pore. One can observe that the counterion current, /+, is always larger in 
the presence of DNA than in the free-DNA case due to the surface induced enrichment and 
that the difference increases at low molarities where the surface conduction mechanism is 
more evident. On the contrary, the coion current, is sensibly higher for the free-DNA 
case at large molarities where the steric hindrance overwhelms the Coulombic repulsion, but 
the difference disappears as the Debye length increases. 

Variations in conductance associated with the presence of DNA are best appreciated by 
considering the relative conductance deviation dehned as the ratio: 

AG = ^^NA-Gfree 

G free 

where Gfree is the conductance of the free pore, while Gdna is the same quantity in the 
presence of DNA. 

Fig. [^illustrates the behavior of such an observable as the electrolyte concentration is varied. 
When the molarity is below a certain threshold value. Me, one observes a positive value of 
AG/G, which corresponds to an enhanced conductance with respect to the free-DNA case. 
Above Me, instead, the presence of the DNA intruder partially hinders the passage of ions 
and determines a reduction of the conductance with respect to its open pore value. Following 
the name convention the two regimes are named enhancement and blockage, respectively. 
We have also tested the importance of the geometry by comparing the relative conductance 
of a straight cylindrical pore with that relative to the double funnel-cylinder system. The 
effect of the inhomogeneity appears to be relevant only at low concentrations (< 0.2M) while 
the critical concentration. Me, results nearly independent on the channel geometry. Fig. 
seems to indicate that the straight pore geometry is more sensitive to the presence of DNA 
since the relative conductance varies faster with decreasing molarity. We ascribe this feature 
to the fact that the absolute value of its conductance is lower than the one associated with 
the double funnel, while the variations due to the DNA intrusion are comparable in the two 
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situations. In Fig. [^we also reported the results obtained using the LEN approximation, 
where one can appreciate that at large molarities the agreement between the data obtained 
in the present work and the LEN is fairly good, whereas at low molarities corresponding to 
larger Debye lengths the LEN theory tends to underestimate the difference in conductance 
between the DNA-free and DNA case. 

The importance of the radius of the pore is stressed in Fig. where the relative 
conductance is shown for Ri = 10, 20, for a reduction of 38% of the nominal charge DNA 
density, while the remaining parameters were the same as in Fig. One can see not 
only that the larger the radius the smaller the crossover concentration, but also that the 
sensitivity decreases with increasing size as one can see from the fact that for a hxed value of 
the molarity the relative conductance of the corresponding to Ri = 20 nm is in general lower 
than the conductance relative to Ri = 10 nm. For the sake of comparison, we also display 
the corresponding results of the LEN approximation and remark that for the crossover value 
appears to be underestimated by the LEN, which predicts that the blockage extends to lower 
molarities with respect to the PNP. 

Fig. [^illustrates numerical results for AG for a non uniform channel of radius Ri = 5nm 
versus salt concentration, varying in the interval 0.05 M a 1 M, corresponding to Debye 
lengths ranging from 0.3 nm to Inm, for different choices of the value of the DNA charge 
(that is for a reduction of 0%, 38% e 70% of the nominal charge DNA —0.61 e/nm?) : as the 
DNA charge decreases the crossover point moves towards lower concentrations in agreement 
with the fact that the conductance is dominated by the surface conduction mechanism 
according to which a larger T^dna determines an enrichment of the counterions and thus a 
larger conductance. In this case the LEN approximation also appears to work quite well as 
compared to the PNP method. 

An interesting aspect regards the role played by the roughness of the surface of the 
DNA molecule. So far, we represented the DNA as a uniform cylinder of radius Rdna and 
neglected its double helix structure. We, now, consider DNA as having a corrugated shape 
of cylindrical symmetry and radius varying along the symmetry axis according to the law: 

27r 

Rdna{x) = rQ +Asm{—x) (17) 

where tq = 1.1 nm , A = 0.5 nm and P = 3.4 nm. Consistently, eqs. (j^ and (j^ are 
modihed in order to take into account the variation of Rdna with position. Fig. shows 
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that the effect of the corrugation is appreciable only in the low concentration regime where 
the surface conduction dominates and is influenced by and the currents are more 

sensitive to the geometrical details of the charge distribution. The same figure also shows 
the comparison between the theoretical results and the corresponding experimental result 
of Dekker and coworkers |3D|. 

A global picture of the conducting properties of DNA-pore system is provided by the 
two dimensional conductance diagram, whose axes are the molar concentration and the 
channel diameter as shown in Fig. For each pore diameter, one determines the critical 
concentration where the relative conductance changes sign, as for instance shown in 
Fig. 11 The plane shown in Fig. is divided in two regions: above the line the ionic 
electric current is reduced since in wider pores the prevailing mechanism of conduction 
is Ohmic and the geometric effect beats the surface conduction; below the line, instead, 
the extra charges made available to conduction within the pore by the presence of fixed 
charges gives rise to the surface conduction and thus to an enhanced current. Each point 
displayed in Fig. [^represents the value of the critical concentration where the crossover 
enhancement/blockage occurs. In Fig [^ for the sake of comparison we have included the 
effect on the variation of the relative conductance of the convective term using the value 


provided by formula (12). The three different lines refer to: a) the channel shown in Fig. 
g b) the same set up but with a corrugated DNA intruder as described by ([^ and c) the 
effect of the convection velocity, Vconv, given by formula (12) on the phase diagram. Fig. 

also shows the effect of DNA convection, which is negligible for radii larger than 7.5 nm. 
Indeed, one can observe that the inclusion of convection determines only a small shift of 
the ’’coexistence line”. On the other hand. Fig. also indicates that the conductance 
diagram is shifted in the direction of higher molar concentrations when one includes the 
possibility of DNA corrugation. In addition. Fig. shows that by increasing the pore radius 
the critical concentration decreases slightly and tends to be independent on the pore size. At 
high concentration the enhancement current can not be recognized. However, it is possible 
to detect the passage of DNA in the pore by the current blockage. 
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IV. CONCLUSIONS 


We have employed a simple one-dimensional representation of a pore-electrolyte-DNA 
system widely used in experiments on DNA translocation. Our focus has been in under¬ 
standing how the ionic currents are altered by the presence of DNA in the pore, modelled as 
a charged cylindrical intruder. The analysis is based on a one-dimensional reduction of the 
three dimensional PNP model which treats the ions in the continuum and the electrostatic 
interactions in a mean held fashion. The resulting effective equations are of diffusive type 
and display the presence of the so-called entropic term, stemming from the variations of the 
geometry of the channel along the direction of the axis, and of a driving electric held due to 
the charges. 

We have found instructive to compare the results of one dimensional PNP model with the 
corresponding results of the TEN theory. The latter in spite of its simplicity can give a hrst 
hint of the behavior of a complex system such as a pore. The LEN calculation shows that 
the crossover from blockage to enhancement is qualitatively reproduced, but quantitatively 
its predictions become more inaccurate especially at low molarities. It is evident that the 
LEN yields crossover values Me for pores of large radius which are too low. We attribute 
this feature, to the fact that the LEN underestimates the effect of the surface charges and 
thus determines a lower conductance. In the future work the results of this one dimensional 
model will be tested against full three dimensional Lattice Boltzmann simulations pTl IB2] . 
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FIG. 2; Comparison between the total electric currents in the case of the double funnel geometry 
(squares) and a straight cylindrical geometry (circles). Inset: Individual ionic currents for the 
double funnel geometry. Circles refer to the case where DNA is present, while squares to the free 
case. Data refer to a channel radius Ri = 5nm, length L = 34: nm and potential difference 0.2 V. 
The wall and DNA surface charge densities are Sq = —0.375 e/nm^ and Sda^a = —0.38175 e/nm^, 
respectively and no convection Vconv = 0. 
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FIG. 3: Relative conductance in the presence of DNA intruder. Comparison between the relative 
conductance versus molar concentration of a cylindrical channel and the 5 stages system of fig. 
corresponding to a radius of the channel i?i = 5 nm and without convection {vconv = 0). The wall 
and DNA surface charge densities are the same as in Fig. [^ . The continuous lines represent the 
corresponding results obtained within the TEN approximation. In the case of the uniform geometry 
the TEN gives a lower value of AG, whereas for the system B the same quantity is larger. 
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FIG. 4: Relative conductances versus bulk ion concentration for channel radii Ri between 5nm e 
20 nm computed within the PNP one dimensional equation and without convection {vconv = 0)- 
The surface charge density of channel wall is —0.375 e/nm^ and T^dna = —0.38175 e/nm^. The 
continuous lines report the corresponding LEN approximation results. 
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FIG. 5: Relative conductance against bulk ion concentration. The channel radius is 5nm , Sq = 
—0.375 e/nm^ and the surface charge concentration density of DNA varies between —0.61 e/nm^ 
and —0.18 e/nm^ and Vconv = 0. The continuous lines report the results obtained via the TEN 
approximation. 
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FIG. 6: Relative conductance versus bulk ion concentration for a pore of radius Ri = 5nm. The 
values of the surface charges are the same as in Fig[^ Circles indicate the PNP results relative to 
the cylindrical DNA, squares those corresponding to the model with corrugations, while diamonds 
are the experimental data of Smeets et al m- 
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FIG. 7: Dependence of the critical concentration on the nanochannel diameter and bulk ion concen¬ 
tration. The surface charge density of nanochannel is —0.375 e/nm^ and = —0.38175 e/nm^ 

. Three cases are presented: non uniform channel (circle), DNA convection (square), corrugated 
DNA (diamond). The convective term Vconv has been obtained by using eq. (12). In the corru¬ 
gated case Rdna{x) was given by eq. 0 with ro = 1.1 nm , A = 0.5 nm and P = 3.4:nm. 
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